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. . . We report the continuum limits of the speed of sound, Ca, and the specific heat at constant 

lO ' volume, Cv, in quenched QCD at temperatures of 2Tc and STc. At these temperatures, Cs is within 

I 2-cr of the ideal gas value, whereas Cv differs significantly from the ideal gas. However, both are 

■ compatible with results expected in a theory with conformal symmetry. We investigate effective 

measures of conformal symmetry in the high temperature phase of QCD. 

, PACS numbers: 12.38.Aw, 11.15.Ha, 05.70.Fh 

in 



o 

(N 



X 



I. INTRODUCTION 
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CO ' It is now a well-established fact that the pressure, P, and the energy density, e, deviate ,1.] in the high temperature 
phase of QCD by about 25% from their ideal gas values at a temperature of about 3Tc, where Tc is the transition 
temperature. Early expectations that e would count the number of degrees of freedom in the QCD plasma through 
.ii— j. I the Stefan-Boltzmann law are belied by the fact that perturbation theory has had great difficulty in reproducing these 
' lattice results There have been many suggestions for the physics implied by the lattice data, the inclusion of 
^ various quasi-particles the necessity of large resummations 10], and effective models |0, Interestingly, there 

■ has been a suggestion that conformal field theory comes closer to the lattice result Q . This assumes more significance 
I ' in view of the fact that a bound on the ratio of the shear viscosity and the entropy density, s, conjectured from 
Q-i, the AdS/CFT correspondence 3 lies close to that inferred from analysis of RHIC data Q and its direct lattice 
^ ■ measurement Q as well as the lattice results of a different transport coefficient 0- 

In this paper we go beyond the measurement of the equation of state (EOS) to thermodynamic fluctuation measures. 
In the pure gluon gas there is only one fluctuation measure, the specific heat at constant volume, Cv Related to this 
is a kinetic variable, the speed of sound, Cg. We report measurements of both of these in the high temperature phase 
of the pure glue plasma, through a continuum extrapolation of results obtained with successively finer lattices. Not 
only do these quantities provide further tests of all the models and expansions which try to explain the lattice data 
on the EOS, they also have direct physical relevance to ongoing experiments at the RHIC in Brookhaven. 

The speed of sound governs the evolution of the fire-ball produced in the heavy-ion collision and hence plays a 
crucial role in the hydrodynamic study of the signatures of QGP. Elliptic flow is one of the most important quantity 
that has been suggested for the signature of QGP formation in the heavy-ion collision experiments. It has been shown 
milllii, 15] that elliptic flow is sensitive to the value of Cs ■ 

In recent years, event-by-event fluctuation of quantities has been of immense interest as signatures of quark-hadron 
phase transition. In order to use fluctuations as a probe of the plasma phase, one has to identify observables whose 
fluctuations survive the freeze-out of the fireball. The evolution of these fiuctuations is sensitive to the values of Cs , 
as shown in |l6j] for the case of net baryon number fluctuation. In [T^ it has been claimed that, within the regime of 
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thermodynamics. 



the ratio of the event-by-event fluctuations of entropy and energy is given by 

{5EY/E^ " {l + Clf 



Re — 



(1) 



and hence provides an estimate of the speed of sound, if i?e turns out to be measurable in heavy-ion collisions. 

The specific heat is, of course, directly a measure of fluctuations. It was suggested in that event-by-event 
temperature fluctuation in the heavy- ion collision experiments be used to measure Cv ■ Also it has been shown in 
that Cv is directly related to the event-by-event transverse momentum [pt) fluctuations. 

The measurement of Cv and Cs also directly test the relevance of conformal symmetry to finite temperature QCD. 
It is common knowledge that QCD generates a scale, Aqcd, microscopically, and thus breaks conformal invariance. 
The strength of the breaking of this symmetry at any scale is parametrized by the /3-function. A recent suggestion is 
that an effective theory which reproduces the results of thermal QCD at long-distance scales may somehow be close 
to a conformal theory. The result of 's^ for the entropy density, s, in a Yang-Mills theory with four supersymmetry 
charges (TV = 4 SYM) and large number of colours, Nc, at strong coupling, is 



-^/Cg^TVj, where /(x) = ^ + -Jc(3)(2.. 



.-3/2N 



So 



32 



and 



(2) 



where g is the Yang-Mills coupling. For the Nc ~ 3 case at hand, the well-known result for the ideal gas, sq — 
4:{Nc — l)7r^r'^/45 takes into account, through the factor — 1, the relatively important difference between a 
SU{Nc) and an U{Nc) theory. 

Of course, at finite temperature there is a scale, T, which appears in, for example, e as a factor of T^. However, 
the strength of the breaking of conformal symmetry must be measured as always, through the trace of the stress- 
tensor. After subtracting the ultraviolet divergent (T = 0) pieces, this is given by the so-called interaction measure, 
A = e — 3P. Thus the ratio, C — A/e, which we call the conformal measure, parametrizes the departure from conformal 
invariance at the long-distance scale '^sij. We shall present data for this quantity in this paper. Note that when A = 
one has C^ = 1/3 and Cv/T^ = 4e/T'*. In addition, if the gas is ideal, then one obtains e/T'* = {N^ - \)tP- jXh. 

The paper is organized as follows. In the next section we present the formalism and lead up to the measurement of 
Cv and Cl on the lattice in Section FlI Bl Some details are given in the appendices. In Section ITTll we give details of 
our simulations and our results. Finally, in Section llVl we present a discussion of the results. Those who are interested 
only in the results may read the initial part of Section ITU and then jump to Section llVI 



II. FORMALISM 



For any theory on the lattice one may compute derivatives of the partition function, Z{y^ T), where V is the volume 
and T the temperature. In particular the energy density, e, and the pressure, P, are given by the first derivatives of 
In^, 



r2 d\nZ{V,T) 



dT 



and 



P = T 



dlTiZ{V,T) 



dV 



(3) 



The second derivatives are measures of fluctuations. For a relativistic gas, where particles may be created and 
destroyed, there is only one second derivative in the absence of chemical potentials, the specific heat at constant 
volume — 



C 



(4) 



A quantity such as the compressibility is not defined in the absence of a chemical potential, since a gas may adapt to 
a change in volume without change in pressure by simply creating or destroying particles. Nevertheless, the speed of 
sound is a well defined quantity. One can see this by using thermodynamic identities to recast the definition in the 
form 



C^ 



dP 

a7 



dP 




^ 5e 




s~df 


J 







where we have used the thermodynamic identity 



dP 
df 



dS_ 
dV 



and — — 



dS_ 
dV 



c^Jt'^ 



e + P 



(5) 



(6) 
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in conjunction with the definition of the entropy density, s, above. Although this work deals exclusively with the pure 
gauge theory, this seems to be a good place to remark that in full QCD without quark chemical potentials all these 
relations go through. However, in the presence of a chemical potential, compressibility is a well defined concept, and 
the speed of sound may be affected by this new physics. We shall deal with this elsewhere. 



A. Energy density and pressure 

In order to set up our notation, and for the sake of completeness, we write down some results which have been 
known since the beginning of lattice thermodynamics. For the pure gauge SU{Nc) theory with the Wilson action, 
the partition function is 



Z{V, T)^ J VUe~^^^\ where S[U] = 2Nc [K^Ps + KrPr] 



(7) 



periodic boundary conditions are assumed, Ps denotes the sum of spatial plaquettes over all lattice sites, and is 
the corresponding sum of mixed space-time plaquettes [s^- If the lattice spacings, {i = s,t) along spatial (s) and 
temporal (t) directions, are related by the anisotropy parameter ^ = a^/ar, then the coupling may be written as 
Kg = 1/^5?! and Kj. = ^/g^. £, = 1 corresponds to the usual case of symmetric lattice spacings, i.e., = a^-- In the 
weak coupling limit, g~'^ 's can be expanded '24'| around their symmetric lattice value g~^{a), 

9^'(as,0 = g-'ia) + c.(C) + 0[g'(a)], (8) 

with the condition Ci{£, = 1) = 0. 

The computation of the first derivatives of the couplings requires b{as) = asdg~^ /dag. With the usual definition of 
the /3-function- 



T3f \ ^ 



— 77: -ai 

127r 



where 



as = .g^/47r. 



one finds b{as) — B{as)/2Tral. Then, 



da 



2TTa^sC 



e 



and 



da. 



27ra2 



K = g-^+^c',, 



(9) 



(10) 



where primes denote derivative with respect to ^. The quantities and c'^ have been computed to one-loop order in 
the weak coupling hmit for SU{Nc) gauge theories p^ . 

If the number of lattice sites in the two directions are Nr and Ns, then T — l/NrUr and V — {Nstts)^- The partial 
derivatives with respect to T and V can be written in terms of the two lattice parameters ^ and Ug , keeping and 
Nt fixed, 



T 



d_ 

dT 



d_ 



and 



3V 



d_ 

dV 



_d_ 

da. 



d_ 



(11) 



Using these derivatives, one obtains [2C 



4 



dKs 

di 



Dg 



dKr 



Dt 



and a^A = (e - 3P) 6iVcCas 



dKg 
da. 



dKr 
da. 



(12) 



A subtraction of the corresponding vacuum (T — 0) quantities lead to Di — {Pi) — (Pq) above, where (Pq) is the 
average plaquette value at T = 0, and the bar denotes division by 3NgNr. The T = quantities are evaluated by 
simulations on symmetric lattices with N^- — Ng ^ oo. In the isotropic limit, ^ = 1, one obtains 



2^4 



Dg - Dt 



{c'gDs + c'rDt 



and 



A 

y4 



Bjas) 
27ra2 



Dg + Dt 



(13) 



Note that A contains B{ag) as a factor, but this explicit breaking of conformal symmetry may be compensated by the 
vanishing of the factor Dg + Dt- To determine the coupling g^, throughout this work, we use the method suggested 
in , where the one-loop order renormalized couplings have been evaluated by using 1^-scheme p7j and taking care 
of the scaling violations due to finite lattice spacing errors using the method in [28| . 
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FIG. 1: di./T'^ as a function of the bare coupling /3 using a non-perturbative (squares) and one-loop order perturbative 
(pentagons) /3-function, B{as). The results agree for /3 > 6.5. The data of Q for A'"t = 8 has been used in this analysis. 



B. The specific heat and speed of sound 



Direct application of the derivatives in eq. Hll|l to the expression for e in eq. Hl'2|l would give an incorrect result for 
Cv, i.e., Cv /T^ would not give the correct ideal gas value {Ae/T'^) in the limit 5 ^ 0. This is because we have chosen 
to work with the variables ^ and as , whereby the dimensions of both T and V come from powers of . Application 
of the derivative formula therefore see false scalings of these quantities. We solve this problem by choosing to work 
in terms of a dimensionless ratio so that the scaling is automatically taken care of. We define 



e 



and 



dT 



(14) 



where C is the conformal measure already discussed. Then, using eqs. (jS) I14|l one can proceed straightaway to write 







' s F e " 


r3 




T3 ^ 3 T4 



and 



C 



3s /T^ 



(15) 



In order to complete these expression, one needs to express F in terms of quantities computable on the lattice. 
We lighten the subsequent formulae by defining the two functions 



Fit as) 



Since C 



F/G, one finds that 



dKs 
da. 



Ds 



dKr 
da. 



Dt 



and 



QN4 



Ds 



dKr 



F = - 



T dF 
G df 



TF dG 
G2" df 



(16) 



(17) 



The derivatives of F and G involve the variances and covariances of the spatial and temporal plaquettes. The complete 
expressions are given in Appendix ^ We also need to write down the second derivatives of the couplings — 



' da. 



B{as) 



K'l 



2.9; 



24 



dK'r ^ Bjas) 
da. l-KcP: 



(18) 
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FIG. 2: Stability of Ds against the statistics in T > (first panel) and T = (second panel) simulations. The coupling 
corresponds to 3Tc on a 12 x 38"^ lattice and the corresponding T — simulation was performed on a 38'* lattice. In both the 
figures we have also plotted the 1-cr error band of the final errors. 



Details of the evaluation of the second order Karsch coefficients c" ^ are given in Appendix El For Nc = 3 their 
numerical values_are c" = -0.298193 and c'^ = 0.333670. 

It is known ^29] that in the weak-coupling limit, i.e., for 5 — s- 0, the dominant contribution to all plaquettes varies 
as g^. So the Z^i's also vary as g^. Hence in this limit A/T"* cx g^ and can be neglected in comparison with e/T^ 
which is {N^ - l)n'^/15 + 0{g^). As shown in Appendix E| T -> in this limit. As a result, C^/T^ 4e/r'* and 
Cg — > 1/3. Note that in any conformal invariant theory in rf + 1 dimensions one has e = c?P, i.e., C — T — 0, and 
hence, by eq. (O, = 1/d and Cy/T^ = {(1+ l)e/T'^. 



C. On the method 



The formalism outlined in Section III Al corresponds to the so-called "differential" method. On the other hand, 
about a decade back a new method called the "integral" method was employed Q to find the equation of state of 
QCD matter. This consists of using the expression for A/T'* in eq. (|13|) while eschewing the expression for e/T* in 
favour of a totally different method for determining P/T^. This proceeds from the observation that P — {T/V)lnZ 
for a homogeneous system. Since lattice computations do not determine the absolute normalization of Z, use of this 
formula requires fixing an additive constant. In this constant was fixed by setting P = in the low temperature 
phase of QCD, slightly below Tc- 

This leads to a conceptual problem. Since the pure gauge phase transition in QCD is of first order, the system is 
not homogeneous at and the method is not strictly applicable there. One makes an unknown systematic error in 
integrating through Tc. This is in addition to a small systematic error due to setting P = just below Tc and the 
numerical integration errors. However, a decade ago, these errors were deemed to be smaller than the scale errors in 
the "differential" method which could give rise even to negative pressure in computations with coarse lattice spacings, 
a — 1/TNr, which were the only ones possible then. 

The two methods must agree in the continuum limit. A reanalysis of the data of showed that the two methods 
agreed to within 5% (i.e., within statistical errors) for T > 2Tc j22| already for Nr = 8, i.e., for a = l/16Tc. In fact a 
criterion for the agreement is straightforward, and follows from the fact the expression for A/T^ is common to the two 
methods. Since the integral method does not need the Karsch coefficients, it allows one to use a non-perturbatively 
determined /3-function B{as). On the other hand, the differential method requires, for internal consistency, that the 
Karsch coefficients and B{as) be obtained at the same order, i.e., at one-loop order in the present state of the art. 

Thus, a comparison between the values of A/T^ extracted for a given N^- using the two techniques would reveal at 
what T the two methods become identical. Then, using asymptotic scaling, one could also give the minimum value 
of Nr which would be required for the same level of agreement as a function of T. Such a comparison is shown in 
Figure Hand demonstrates that a bare coupling of /? > 6.55 already suffices. 

The use of the "differential" method does require small lattice spacings, but this is currently quite feasible with 
modern computers. The advantage is that the formalism of Section [ll BI can be used to compute second derivatives of 
the free energy as operator expectation values. If the EOS were to be evaluated by the "integral" method, then these 
quantities could only be evaluated by numerical differentiation 23] , which is a less attractive solution. 
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FIG. 3: Dependence of Ds against the spatial size, Ns, of the T > lattice for 2Tc and 22^ lattice for the T — computation 
(left panel). We have shown a fit to a constant through the three largest lattices. Dependence of Ds on Ns for the T = lattice 
when the 2Tc computation is performed on an 8 x 18^ lattice (right panel). A fit to a constant with the two largest sizes is also 
shown. 
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FIG. 4: Dependence of e/T* (squares) and SP/T* (pentagons) on for T = 2Tc (left panel) and T = 3Tc (right panel). 

The 1-a error band of the continuum values has been indicated by arrows (for e/T*) and lines (for 3P/T"^). 
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TABLE I: The coupling (/3), lattice sizes (Nr x A'^^), statistics and symmetric lattice sizes (A*) are given for each temperature. 
Statistics means number of sweeps used for measurement of plaquettes after discarding for thermalization. 



Nr 


?,QNlN^Vos,s 






10 


(0.2 ± 1.0) X 10"^ 


(0.2 ± 1.0) X 10"^ 


(0.1 ± 1.0) X 


10"^ 


12 


(0.7 ±4.1) X 10"* 


(0.7 ±4.1) X 10~* 


(0.4 ±3.5) X 


10"* 


14 


(0.1 ± 1.4) X 10"^ 


(0.1 ± 1.5) X 10"^ 


(0.1 ± 1.3) X 


10"^ 



TABLE IL Contributions of different covariances, with their respective errors, are tabulated for different N-r at 2Tc. This shows 
all the contributions are totally negligible. 
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FIG. 6: In the left panel we show the continuum values of 4e/T* (pentagons) and Cv /T'^ (squares) against T /Tc. In the right 
panel we show the temperature dependence of the continuum extrapolated values of . 
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FIG. 7: Continuum values of the conformal measure, C, (squares) and the one-loop order /3-function (circles) evaluated at the 
scale 2ttT. For the T = gluon condensate see the text. 



III. SIMULATIONS AND RESULTS 



The simulations have been performed using the Cabbibo-Marinari pseudo-heatbath algorithm with Kennedy- 
Pendleton updating of three SU{2) subgroups on each sweep. Plaquettes were measured on each sweep. For each 
simulation we discarded around 5000 initial sweeps for thermalization. We found that integrated autocorrelation time 
for the plaquettes never exceeded 3 sweeps. In Tabled we give the details of our runs. All errors were calculated 
by the jack-knife method, where the length of each deleted block was chosen to be at least six times the maximum 
integrated autocorrelation time of all the simulations used for that calculation. 

Of all the quantities which go into determining the equation of state, namely Dg, Dt and Dg — Dt, we found that Dg 
was the smallest. Hence control over the errors of Dg was the most stringent requirement on the amount of statistics 
needed. In Figure [21 we show the stability of Dg against the finite temperature and zero temperature statistics for 
the case where we have minimum statistics, namely for the 3Tc run on the 12 x 38^ lattice and the run at the same 
coupling on the 38"^ lattice. Note that the plaquette values are of the order of unity, and the first four digits cancel in 
computing Dg. Thus, the control over errors shown in Figure El was due to our reducing the errors in the plaquette 
variables to a few parts per million. 

We checked the stability of Dg against the spatial size of the T > lattice. The left panel of Figure (31 displays the 
dependence of Dg on Ng for Nr — 4 when the temperature is 2Tc and the symmetric lattice size is 22*. We have also 
shown a fit to a constant using the data on the three largest lattices. These fits pass through the data collected on the 
Ng = {T/Tc)Nt + 2 lattice. We have checked that this condition holds for 3Tc also. This observation is consistent with 
the results of earlier investigation of finite size effects for T > 30] and motivated our choice of Ng = {T/Tc)Nr + 2. 

We also investigated the dependence of Dg on the size of the symmetric Ng lattice used for the T — subtraction. 
The right panel of Figure [31 exhibits the dependence of Dg on Ng for a run at 2Tc on a 8 x 18'^ lattice. Dg is seen to 
be constant for Ng > 22. In view of this we have used 22* as our minimum T = lattice size, and scaled this up with 
changes in the lattice spacing. 

We evaluated the contributions of the terms containing different covariances of the plaquettes and found them to 
be totally negligible, as shown in Table ITll It is worth noting that a previous computation in SU(2) pure gauge theory 
close to Tc found significantly larger, and clearly non- vanishing, variances of the plaquettes (31|. This suggests that 
the variance terms might make significant contributions to the specific heat and speed of sound closer to the softest 
point of the EOS. 

Following the observation illustrated in Figure ^ all our continuum extrapolations have been done with lattice 
spacings which are at least as small as that at /3 = 6.55. For T ~ 2Tc this leaves three values of Nr from which an 
continuum extrapolation linear in oc l/N^ can be performed. However, at T = 3Tc, the continuum extrapolation 
has been performed with two values of Nr- This completely fixes the two parameters of the linear extrapolation 
to the continuum and the error in the continuum extrapolated value is by definition zero. In this case the error in 
the continuum value was estimated by three methods. Two of these consisted of first making the best fit using two 
parameters, and then keeping one fixed while allowing the other to vary in order to make an estimate of the error in 
that parameter. We also made extrapolations using the upper end of one error bar and the lower end of the other. 
This last procedure gave the maximum errors in the continuum extrapolated values, and we choose to quote this, 
since it is the most conservative error estimate. 
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T/Tc 






p/t" 




Cv/T'-'' 




2.0 


8.1 (1) 


3.3 (3) 


1.0 (1) 


4.2 (3) 


15 (1) 


0.30 (1) 


3.0 


7.0 (1) 


4.2 (5) 


1.3 (2) 


5.6 (6) 


18 (2) 


0.30 (2) 



TABLE III: Continuum values of some quantities at the two temperatures we have explored. The numbers in brackets are the 
error on the least significant digit. The value of the t'Hooft coupling NcQ^ is computed at the scale 2txT using the Tc/Aj^ 
quoted in [2^. 



In Figure^ we present our results for e and P. At STc the continuum limit values of e/T^ and P/T'^ differ from 
their respective ideal gas values (87r^/15 and 87r^/45 respectively) by about 24%, and by about 40% at 2Tc. This is 
consistent with previous measurements of the equation of state at these temperatures. At 2Tc our results differ from 
the ideal gas value by almost 1-a. The actual results are also collected in Table Hill 

Similar continuum extrapolation of Cv/T^ and are shown in Figures and Figures El respectively. Continuum 
extrapolated values of different quantities for two temperatures are tabulated in Table Hill It can be seen that although 
e/T"* and P/T^ differ significantly from their ideal gas values, is quite close to 1/3, being within 2-3 a. Similarly, 
it can be seen that Cv/T^ is completely compatible with Ae/T'^, a fact also illustrated in Figure [7| However, its 
deviation from the ideal gas value is seen to be more significant. 

The generation of a scale and the consequent breaking of conformal invariance at short distances in QCD is quantified 
by the /3- function of QCD, and at long distance in the finite temperature plasma by the conformal measure C. In 
Figure El we plot these two together. The T = gluon condensate, (-F^}o gives a measure of the non-perturbative 
strength of the breaking of conformal invariance at T = 0. A recent estimate 33] gives {F^)q/A±^ = 10.5 ±0.5, where 
A-jjs is compatible with the value used to set the temperature scale. It is clear that conformal invariance is better 
respected in the finite temperature effective long-distance theory than at the microscopic scale. 

An alternative measure of conformal symmetry breaking for T > could be the dimensionless quantity 

Clearly its use is limited to making the trivial point that the scale T breaks conformal symmetry explicitly, reflected 
by the growth of C This misses physics on both the weak-coupling and the strong coupling AdS/CFT sides. On 
one side is the likelihood that the weak coupling expansion, which is conformal, becomes more accurate at sufficiently 
high temperature — although A increases with T, A/T** becomes a small part of e/T^, P/T^ or s/T'^. On the other 
hand is the possibility that AdS/CFT is accurate because the T dependence can be taken into the dimension of 
operators, and the CFT shows up in non-trivial dimensionless quantities such as the function / in eq. |2Jl. C allows 
us to test the accuracy of these possibilities, thereby capturing the interesting facts shown in Figure |B1 whereas C 
misses them. 



IV. DISCUSSION 



In this paper we have determined the continuum limits of the specific heat at constant volume, Cv, and the speed 
of sound, Cg, in the continuum limit of the pure gluon plasma. In the process we have also recomputed the EOS, 
i.e., the pressure, P, and the energy density, e, by a method which has not been used earlier to obtain the continuum 
limit. Our results are collected together in Table HTTl It is clear from this that the plasma phase is not far from the 
conformal symmetric limit, in which — 1/3 and Cv/T^ — Ae/T^. This is also illustrated in Figures |S| and |7| The 
latter figure further shows that conformal symmetry is more nearly realized on length scale of order T than one would 
expect either from the /3-function of QCD or from the gluon condensate (P^)o- 

We have extended the operator method for the EOS in order to obtain these results. We have limited our com- 
putations to lattice spacings where the method is guaranteed to work because of the observed scaling of results with 
the correct QCD beta function, as shown in Figure ^ The main advantage of this method is the relatively noise- free 
determination of Cg and Cy 

There is as yet no unambiguous determination of these quantities from experiments. Hydrodynamic analyses of 
RHIC data routinely use = 1/3 in the plasma phase; our results show that this receives a correction of up to 10% 
in QCD. It would be interesting to check the sensitivity of the RHIC data to the variation in and compare the 
bounds there with our lattice results. Little is known yet about Cv The only estimate that we are aware of M uses 
the data on fluctuations from ^ to obtain C^/T^ = 60 ± 100 for T = 180 MeV. 
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FIG. 8: The equation of state of QCD matter. The diagonal line denotes possible EOS for theories with conformal symmetry. 
The circle on the diagonal denotes the ideal gluon gas, whose EOS in this form is temperature independent. The ellipses denote 
66% error bounds on the measured EOS. The ratio of the axes is a measure of the covariance in the measurements of e/T* 
and P/T*, which is about 90%. The wedges piercing these ellipses have average slope Cf, and the opening half-angle of these 
wedges indicate the error in Cj. 



Our data on the entropy density, s/T^, shown in Table ITTll can be used to test the quantitative predictions of the 
strong-coupling expansion of the A/" = 4 supersymmetric Yang-Mills theory in (3J . Eq, 10) should interpolate the ratio 
s/sq between the strong-coupling limit of 3/4 and the weak-coupling limit of 1. It is clear that at T = 2Tc this ratio 
is below 3/4, and hence cannot be described by the computation in 0]. At 3Tc, the ratio is above 3/4 and there 
is quantitative agreement with the formula of Q using the value of the t'Hooft coupling quoted in Table Hill It is 
interesting that this particular strong-coupling theory fails at temperatures where weak coupling theories also break 
down. 

A partial summary of our results is illustrated by plotting the equation of state as P/T* against e/T'*, as in Figure 
|S1 In this plot, the ideal gas for fixed number of colours is represented by a single point which is independent of 
T, and theories with conformal symmetry by the line e — 3P. Pure gauge QCD lies close to the conformal line 
at high temperature, as shown. One expects the EOS to drop well below this line near Tc, since the theory then 
contains massive hadrons (glueballs in pure gluon QCD) with masses in excess of Tc. Thus one expects the origin to 
be approached almost horizontally as T ^ 0. We will present data closer to, and below, Tc in a forthcoming paper. 

SG would like to thank Jean-Paul Blaizot, Kari Eskola, Rajesh Gopakumar, Past Huovinen, Toni Rebhan and Vesa 
Ruuskanen for discussions. RVG would like to thank the Alexander von Humboldt Foundation for financial support 
and the members of the Theoretical Physics Department of Bielefeld university for their kind hospitality. SM would 
like to thank Rajarshi Ray for discussions and for his help in calculating the derivatives of Karsch Coefficients. SM 
would also like to thank TIER Alumni Association for partial financial support. 



APPENDIX A: DERIVATIVES OF F AND G 



From the definition of F in eq. IjlGI) we obtain 

'dK' 



dF 



da. 



(9a,, 



9a, 



dKr 
9a, 



(Al) 
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Since the plaquette operators do not explicitly depend on ^ (or a^), we can easily take the derivative — 

D[ = D,{S') - m - Po)S') = -6N,NrN^ [K'^a^a + K^rA , (A2) 

where ag^i = {DgDi) — {Ds){Di) and Ur^i = (DtDi) — {Dt){Di) are the variances and covariances of the plaquettes. 
Next, using Eq ifTUIl and Eq (fTH|l . the ^ 1 limit becomes 



dF 

Similarly we find 



27ra? 



[A - Ds] - 6N,NrN: 



■ B{a,) 
27ra2 



(A3) 



dG_ 



9^ 



9^ 



.(A4) 



Since Di ~ and (Tij- ~ in the limit g ^ Q, dF/d£, but dG/d£, remains finite. As a result, F ^ 0, in this 
limit, as it should, even when the variances of the plaquettes are taken into account. 



APPENDIX B: SECOND ORDER KARSCH COEFFICIENTS 



Expressions for Cs{£,) and Ct-(^) as given in |23| are 
'N^ - 1 3 



^7(0 + ^F/7V(0 + 0.010245 



Here 6^ = gin xi + sin X2 + sin 2:3, and 



6(C'+ 62)1/2 J^^^(^2^^2)l/2y 



^3(0 - e 



2\" /-^g . 2 .2 e2+2&2 

d a;sm Xi sm X2 



/4(^) = ^ 0^ y d3^sin2 2a;i 



53(^2.^^2)3/2' 
^2 + 2&2 



53(^2+62)3/2 



/6(^) = e'(^|y|d3^6-l(C2+62)-3/2cos2xi 



(Bl) 



(B2) 
(B3) 
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h ha hb h h h h h FIN 


Integrals 


0.750000 0.929600 0.119734 0.103289 0.478934 0.250000 0.206578 0.309867 0.0 


1-st derivatives 


0.440133 0.208546 0.190133 0.033774 0.065779 0.309867 0.238384 0.411188 0.003166 


2-nd derivatives 


-0.518412 -0.663270 0.030921 -0.088777 -0.262319 -0.101321 -0.146561 -0.159106 -0.014471 



TABLE IV: Values of the integrals Ix[C) 's and FIN{^) and their derivatives with respect to ^ at ^ = 1. 



Limits of all the above integrals are [0,7r/2]. And 



'-7r/2 

FIN{^) = DIV{^) - DIV{1) 



(B4) 



The integral DIV{^) is infrared divergent. However, what actually needed here is FIN{^). FIN{^) is not divergent 
as it is constructed by subtracting the integral DIV{1), having the same infrared divergence as that of DIV{^), from 
DIV{£,). The derivatives of FIN{£^) are just the derivatives of DIV{^) which are not divergent. We calculated 
FIN"{1) by direct numerical integration and also by taking a derivative of FIN'{^), at ^ = 1, numerically. We found 
that values obtained from both the methods are consistent. 

The numerical values of all these above integrals and their derivatives with respected to ^, at ^ = 1, are tabulated 
in Table Hvl Using these values we can obtain ^ derivatives of the Karsch coefficients at ^ = 1 



Nl-l. 
32A^2 

327V2 

32N^ 
-N^ - 1 



327V2 



0.586844-1-0.000499 
0.691216- 0.005649 
0.586844 + 0.005306 
1.038595- 0.001044 



(B5) 



The values of all the integrals and their first derivatives as well as the regular Karsch coefficients match with their 
respective values mentioned in [2^ . 
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